function [A] = tensorProduct ( a , b , dim)

A = zeros (dim,dim);

for i=1:dim
    for j=1:dim
        A(i,j) = a(i) * b(j);
    end
end


end
